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The presence of the long-range Coulomb force in dense stellar matter implies that the total charge 
cannot be associated with a chemical potential, even if it is a conserved quantity. As a consequence, 
the analytical properties of the partition sum are modified, changing the order of the phase tran- 
sitions and affecting the possible occurrence of critical behaviours. The peculiar thermodynamic 
properties of this system can be understood introducing a model hamiltonian where each charge is 
independently neutralized by a uniform background of opposite charge. Some consequences on the 
characteristics of mixed-phase structures in neutron star crusts and supernova cores are discussed. 
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I. INTRODUCTION 

Among the major quests in investigating the proper- 
ties of compact-star matter, is the underlying equation 
of state. In particular, exotic phases (the so-called pasta 
phases) are predicted in the outer crust of neutron stars 
as well as in the inner core of exploding supernovae. In 
this latter physical situation matter attains temperatures 
of several MeV and the possible existence of such 

dishomogeneous structures is known to have some rele- 
vance for hydrodynamic properties of stellar matter and 
neutrino transport [4, 5] . 

The interplay between the short-range nuclear force 
and the long-range non-saturating Coulomb interaction 
is responsible for the existence of such phases. First cal- 
culations of the composition of the stellar crust were done 
using the Liquid Drop Model or in the Thomas- Fermi ap- 
proximation 0, H, with an a-priori assumption 
on the different species constituting the matter d, , or 
on the shape of the pasta structures 0, H, : the pre- 
ferred shape was essentially determined by the balance 
between surface and Coulomb energies. In these stud- 
ies, as well as in more recent ones along the same lines 
[ill , m, [H, 13 1 charge neutrality is imposed only glob- 
ally while charge fluctuations are allowed on any (small) 
scale, meaning that the effect of the electron background 
is neglected. Then the problem reduces [l^ to a sim- 
ple application of Gibbs phase-coexistence conditions in 
multi-fluid systems As a consequence, the transition 
from the homogeneous outer core to the pasta phase(s) is 
systematically considered as first order in all these works. 

However, the conditions to define a thermodynamic 
limit in systems involving non-saturating forces are not 
trivial and phase transitions may have specific proper- 
ties analogous to those discussed for finite systems |17| . 
In this paper we show that the suppression of one ther- 
modynamic degree of freedom due to the divergence of 
the Coulomb energy-density for any net charge at the 
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thermodynamic limit implies the disappearance of phase 
transitions or a modification of their order. Specifically, 
we show that the crust-core transition is not first order 
but continuous. 



Another limitation of the first seminal works on the 
phase structure of stellar matter is that most calcula- 
tions were done in the mean-field approximation, which 
is known to be especially poor in describing phase transi- 
tions. Recent calculations with semiclassical models be- 
yond the mean field [1, Q naturally include thermal fluc- 
tuations and allow the rearrangement of the proton dis- 
tribution due to the Coulomb field under the explicit con- 
straint of charge neutrality. The inclusion of such fluctu- 
ations was expected to lead to an increased matter opac- 
ity to neutrino scattering with important consequences 
on the supernova cooling dynamics 5, 18]. Such a co- 
herent neutrino-matter scattering is not only expected 
at low temperature, where the pasta-core transition was 
supposed to be first order, but even more in the case 
of the occurrence of a critical point in the post-bounce 
supernova dynamics with the associated phenomenon of 
critical opalescence [H, [l^. Surprisingly, the expected 
increase in the static form factor is not observed in finite- 
temperature molecular-dynamics calculations [20| . indi- 
cating a fluctuation suppression respect to the simplistic 
scenario of a first-order phase transition. 



On the basis of simple scaling arguments, we show in 
this paper that critical behaviors can survive up to the 
thermodynamic limit only if the particle-density fluctu- 
ations associated with the mixed pasta phases do not 
correspond, close to the critical point, to charge-density 
fluctuations. In turn, this would imply a strong increase 
with temperature of free-proton drip or, alternatively, a 
very strong polarization of the electron field, leading to a 
complete charge screening of the pasta structures. Since 
both effects are likely to be unphysical, critical fluctua- 
tions and the associated critical opalescence are expected 
to be quenched. 
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II. THERMODYNAMIC FEATURES OF 
(PROTO-) NEUTRON STARS 

A. Thermodynamics of charged systems 

Let us consider the dense matter in neutron-star crusts 
and supernova cores formed of electrons (e), neutrons (n) 
and protons (p). The microscopic Hamiltonian reads 

H ^ Hnp+ Ke + Vee + Vpp + Vep , (1) 

where Hnp is the nuclear strong interaction including the 
nucleon kinetic energy, Kg is the electron kinetic term 
and Vii' is the Coulomb interaction between different 
types of particles {i e {e,n,p}) : 

v.. = ^ / '-^P^drir- , (2) 
1 + J |r - P 

where a is the fine-structure constant and Pi{r) is the 
local density of the particle of type i and charge qi = ±1. 

Since the Coulomb field is a long-range interaction, 
the existence of a thermodynamics, i.e. the conver- 
gence of a thermodynamic limit and the equivalence be- 
tween statistical ensembles, is not guaranteed fT7|. Let 
us first consider the canonical ensemble with densities 
Pi — Ni/fl, where is the volume containing a num- 
ber Ni of particles i ~ n, p or e. If the net charge 
density Pc ^ Pp — pe is not strictly zero at the ther- 
modynamic limit, the Coulomb energy per unit volume 
{Vc)/^ oc p^ri^/^ diverges. The strict charge neutral- 
ity demanded by the thermodynamic limit guarantees 
the additivity of thermodynamic potentials as well as en- 
semble equivalence. Indeed, the monopole contribution 
to the Coulomb energy between two separated neutral 
systems is identically zero. The longer-range multipole 
interaction is the dipole-dipole term E (x Di ■ D2/ , 
where D is the dipole moment in each subsystem. This 
moment is proportional to NcR, where Nc is the local- 
ized charge number and R the maximum size of the actual 
dipole. To avoid the divergence of the Coulomb-energy 
density, Nc must scale at maximum proportionally to 
R^, indicating the convergence of the dipole part of the 
inter-system energy density, and the additivity of the two 
subsystems at the thermodynamic limit. This reasoning 
can be extended to all multipole-multipole interactions. 
The Coulomb interaction between neutral systems then 
behaves like short-range interactions [l^, leading to en- 
semble equivalence according to the standard demonstra- 
tion of the Van Hove theorem. 



B. Thermodynamic consequences of charge 
neutrality 

It is important to remark that the strict neutrality con- 
straint discussed above is very different from the trivial 
condition of global charge neutrality of the stellar object. 



which is universally recognized. We have just argued that 
neutrality is a necessary thermodynamic condition: this 
means that charge dishomogeneities can appear only at 
the microscopic level and no locally charged domain can 
exist at a macroscopic scale, even if this scale is small 
compared to the size of the star. From a practical point 
of view, this means that in hydrodynamical calculations 
the condition Pp^f) — Pe{r) has to be imposed at every 
location r, and that this neutrality constraint must be 
explicitly applied at the scale of the Wigner Seitz cell in 
neutron-star crusts. This can have important effects on 
the phase structure of matter even if the electron back- 
ground is assumed to be uniform, as we develop in great 
detail in the next sections. 

The important consequence of the charge-neutrality 
constraint is that the canonical free-energy density / 
is defined only for pc = 0. Hence f (T, pn, Pp, Pe) = 
f {T,pn,p) with p — [pp + Pe)/2 and the chemical po- 
tential associated with pc can not be defined since the 
free energy is not differentiable in the total-charge direc- 
tion. We can notice that this can also be deduced in 
a grand-canonical treatment. Indeed the divergence of 
the Coulomb-energy density for pc 7^ forces the grand- 
canonical partitions to fulfill the constraint pp = pe ex- 
actly for all couples of chemical potentials pe and p,p. If 
the chemical potentials associated with pc and p are in- 
troduced in the form pc = [pp — /ie)/2 and p = pe + IJ-p 
respectively, the grand potential per unit volume g results 
independent of pc, so that g{T, pn, p, Pc) — 9{T, Pn, p)- 
Indeed, densities are grand-potential derivatives : the 
condition on the total charge pc = then corresponds to 
a fiat grand-potential in the chemical-potential direction 

This suppression of one degree of freedom arises from 
the thermodynamic limit and should not be confused 
with an additional constraint as it is often done in the 
literature, where the two conditions of charge neutrality 
(expressed as Up = Ue) and /3 equilibrium (expressed as 
Pp = Pe) are treated on the same footing [15| or con- 
fused [1^. Additional constraints, such as constant par- 
ticle fraction or chemical (/3) equilibrium, may or may 
not be realized in the supernova evolution; they are re- 
strictions of the accessible states and do not affect the 
thermodynamical properties, which are state functions. 
Conversely, charge neutrality has to be fulfilled for each 
(macroscopic) physical state. This changes the the num- 
ber of degrees of freedom of the thermodynamic poten- 
tials, which directly affects the thermodynamics. In par- 
ticular, phase coexistence can occur only between two 
neutral phases. The practical consequences of this ther- 
modynamic requirement are detailed in the next sections. 



C. Quenching of mean- field instabilities in stellar 
matter 

These considerations are especially relevant when 
phase transitions are concerned. Phase transitions should 
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FIG. 1: (Color online) Minimal free-energy curvature of ho- 
mogeneous matter at the fixed finite temperature T = 10 MeV 
in the mean-field approximation, obtained with the Sly230a 
interaction. Without Coulomb interaction : nNM (neutral 
Nuclear Matter). Switching on Coulomb interaction : cNM-|-e 
(charged Nuclear Matter neutralized by an electron gas). Dot- 
ted line : nuclear matter with electrons at T = 0. 



be analyzed by looking at the derivative and curvature 
properties of / in the three-dimensional state-variable 
space (T, p„,p). In particular, a first-order phase tran- 
sition at a given temperature identifies with the linear 
behavior of the free-energy density / between at least 
two points A = {p^,p^) and B = {p^,p^), meaning 
that all intensive parameters in A and B are equal (l6j . 
However, since the chemical potential associated with the 
total charge is not a thermodynamic variable, the condi- 
tion fj,^ = p,^ docs not imply as it is often assumed in 
the literature [H, ll^, ll5|, |2l| that the chemical potentials 
fie and fip are both identical in the two phases {fif = pf 
and Pp = Pp). The difference in the chemical potential 
of charged particles is counterbalanced by the Coulomb 
force: as some electrons move from one phase to the other 
driven by the chemical-potential difference, the Coulomb 
force reacts forbidding a macroscopic charge to appear. 
To illustrate this point, we consider a mean- field approx- 
imation, where the free energy of the system is defined 
as the sum of independent baryon ft and electron fe free 
energies, 

/(T, Pn, p) = fb{T, Pn,Pp^ p)+ fe{T, Pe ^ p) ■ (3) 

To spot a phase transition we can study the convexity 
anomalies of / looking at the curvature matrix 



C 



dpn/dpn dpn/dpp 

dpp/dpn dpp/dpp + dpe/dpe 



(4) 



where we have introduced the chemical potentials pn = 
dfb/dpn , Pp = dfb/dpp and pe = dfe/dpe- The addi- 
tional term = dpe/dpe in the matrix modifies the 
stability conditions with respect to the nuclear-matter 
part, i.e. to the curvature of fb- In general, we can expect 



a quenching of the instability: since the electron suscep- 
tibility is always positive, the instability conditions 
corresponding to Det C < or tr C < are more difficult 
to fulfill. Moreover, due to their very small mass, elec- 
trons are almost always a relativistic degenerate Fermi 
gas: 



9Pe 



(5) 



Since is large at subsaturation densities, the quench- 
ing is expected to be strong. A quantitative application 
is shown in Fig.l, which displays the minimal eigenvalue 
of the free-energy curvature matrix eq. Q as a function 
of the isospin asymmetry ps/pf, = (/?„ - Pp)/{pn + Pp), 
evaluated with the effective Sly230a Skyrmc interaction 

At the thermodynamic limit, the nuclear free-energy 
density appearing in eq.([31) is defined for each tempera- 
ture and each point (/9„, pp = p) in the density plane by 
a Legendre transform of the Gibbs free-energy density 



(6) 



9=n,i 



In the mean-field approximation, this latter can be 
approximated using independent-particle averages only 



g^-T\nz^--[{W)-{H) 



(7) 



In this expression, zq is the grand-canonical partition sum 
of a non-interacting Fermi gas of free protons and neu- 
trons with effective mass m*, and chemical potentials 
/i* shifted of the value of the mean field /i* = pq — Uq. 

{W) = J2q{P^ /'^''^q + Uq) is the average single particle 
energy, (H) = Wi. is the average energy with the chosen 
Skyrme functional, and the local mean field and effective 
mass are defined as partial derivatives of the energy den- 
sity Ti. with respect to the particle densities pq and the 



momentum densities Tq = {p^) /h^ 



m 

dpq 



dH 



2m* dTq 



(8) 



All averages are thermal averages over the modified 
single-particle Fermi distribution 



"'^^^ 1 + exp{Pipy2m*q - p*q)) 



(9) 



From Figure 1 we can see that ordinary neutral nu- 
clear matter (nNM) shows a wide unstable (spinodal) 
region characterized by a direction of negative convexity 
of the free energy. In this region, any density fluctuation 
of homogeneous matter is spontaneously amplified, and 
leads to phase separation. This result is drastically mod- 
ified considering charged nuclear matter neutralized by 
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an electron gas (cNM+c) : the free energy is everywhere 
convex, and this stays true at zero temperature (dashed 
Hne) . Since the extension of the spinodal region mono- 
tonically decreases with increasing temperatures, this lat- 
ter result implies that stellar matter does not present any 
thermodynamic instability in this model. This quench- 
ing is due to the fact that matter dishomogeneities cause 
charge dishomogeneities, which cannot be fully screened 
by the electron gas because of its high incompressibility, 
and are thus forbidden at any (macroscopic) scale. As a 
consequence, the different structures in the stellar crust 
(free nucleons and nuclei, bubbles or pasta) have to be 
considered as intrinsically microscopic objects, and can- 
not be treated as different phases in coexistence, obeying 
standard phase equilibrium rules, as it is often done in the 
literature [13, [HI ■ We will see in section IIIFI that such 
microscopic charge fluctuations naturally appear beyond 
the mean field, and they indeed dominate the density re- 
gion between the liquid core and the gazeous outer crust. 
Such structures are however associated to an important 
energetic cost, which, as we will demonstrate, leads to an 
expansion of the dishomogeneous phase. 

If the qualitative effect of electron quenching of the 
instability is independent of the model, the net quan- 
titative effect on the spinodal zone depends on the pa- 
rameters of the nuclear interaction: the different Skyrme 
forces we have analyzed [23] all lead to a complete in- 
stability suppression, while a small spinodal region still 
subsists at zero temperature in NL3 RMF stellar matter 
calculations p6l |. 

D. Suppression of first-order transitions in stellar 
matter 

In order to expose how the charge-neutrality con- 
straint affects the thermodynamics of star matter, we 
now study the core-crust transition in two different ap- 
proaches. One corresponds to the physical situation, for 
which electro-neutrality has to be ensured by the cor- 
relation between proton and electron density distribu- 
tions. For the other approach, we consider a model sys- 
tem where each charged particle species is independently 
neutralized by a uniform background of opposite charge. 

The modified Hamiltonian reads 

H' = Hnp + k, + VI, + v;^ + v;; , ( lo) 

where the interaction with the opposite background 
charge has been subtracted to the various Coulomb terms 

with 5pi{r) — pi{r) — (pi). This system does not present 
any divergence so that the free energy f'{T,pn,Pp,Pe) 
can always be defined. When pp = p, the introduced 
backgrounds cancel out exactly. In this case, the canoni- 
cal ensemble of our initial stellar problem is equivalent to 



the canonical ensemble of the modified Coulomb system 
along the pp = p, surface : 

/(T, p„, p) = f'{T, p„, Pp = p,pe = p) . (12) 

The two pressures are connected by the relation 

g{T,Pn,p) ^ g'{T,p[^,p'p,p'J , (13) 

where the chemical potentials /i^ — df /dpi of the differ- 
ent species i = n, p, e, are linked to the physical chemical 
potentials by the relations /i„ = /ij^, = /i^ -I- /Xg. Let 
us however notice that this equivalence breaks down in 
phase coexistence, since a phase mixture in the neutral- 
ized model belongs to the physical model only if pc = 
in each phase separately. 

In the neutralized Coulomb system, all chemical poten- 
tials and densities are allowed. This model may present 
phase transitions according to the usual phenomenology 
of multifl-uid systems 16]. In particular, first-order phase 
transitions are characterized by the equality of all inten- 
sive parameters between two points A' and B' of the 
particle-density space: (pf > Pp : pf ) and (pf , pf , pf ). 
The linear behavior of /'(T, pn, Pp, Pe) between A' and B' 
corresponds to a single point in the intensive-parameter 
space (T, p'„, p^, p'g) with different densities (i.e. pressure 
derivatives) on both sides. 

An illustration of this approach is given by figures 2 
and 3, presenting mean- field calculations of such neu- 
tralized model, where nuclear and electron contributions 
are treated independently. The strong hamiltonian is 
parametrized by the Skyrme Sly230a interaction ^2§\ and 
the nuclear-matter thermodynamics is evaluated within 
the mean-field approximation, as in the previous section 
[H, m, 1131 • Until a critical temperature, nuclear mat- 
ter (nNM) presents a first-order liquid-gas phase transi- 
tion. Electrons, conversely, are treated as a single ho- 
mogeneous Fermi gas. This results in a violation of the 
relation pp — p, in each phase, since we have pp > (<)pe 
in the dense (dilute) phase. In other words, pc = Pp — Pe 
is an order parameter of the phase coexistence. This is 
permitted within the neutralized model where proton and 
electron charge are independently cancelled in each phase 
according to eq. (|lip , suppressing the Coulomb contribu- 
tion. 

Figure 2 emphasizes the role of pc as an order param- 
eter for a phase coexistence occurring in the neutralized 
model. In the neutron-star phenomenology, this neutral- 
ized system corresponds to nuclear-matter calculations 
[l^ which do not explicitly impose the neutrality con- 
straint on the scale of the Wigner-Seitz cell. The high- 
density phase represents the homogeneous liquid core, 
while the low-density one is the unbound vapor. Let us 
consider a system such that pc — 0. The corresponding 
path is represented on Figure 3 as a projection in the 
chemical-potential plane (pe,Pp), for different values of 
the neutron chemical potential p„. Starting from a high- 
density homogeneous system, and decreasing the density, 
the liquid-gas phase coexistence is reached at point L. 
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FIG. 2: (Color online) Net charge density pc = —dg'/d^.'c 
as a function of the proton and electron chemical potential 
for a neutralized Coulomb system at fixed temperature and 
neutron chemical potential, with the Sly230a effective interac- 
tion. The solid line gives the physical constraint pc = cte = 0. 
The neutralized model presenting a first order phase transi- 
tion with a pc discontinuity, the physical constraint forces the 
system to follow the coexistence line, leading to a continuous 
transition for the physical system. 
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FIG. 3: Upper left : representation of the constraint pp — pe 
as paths in the chemical-potential plane (^e, p,p). Upper right 
: for each given neutron chemical potential, proton density is 
determined by pp. Lower part : electron density is uniquely 
determined by pe, according to eq.(l5]). /i„ = —17.17 MeV cor- 
responds to the coexistence zone of symmetric matter, higher 
Pn gives a matter increasingly neutron rich. 



Afterwards, the global requirement Pc = has to be ful- 
filled by phase coexistence between two phases for which 
Pc ^ 0. An exemple is given by points A' and B' , which 
obey the phase-equilibrium condition of equal intensive 
parameters. This is a first-order phase transition. Similar 
features are found for any fixed under-critical tempera- 
ture. 



However, this does not correspond to the physical sit- 
uation, where electro-neutrality imposes pc ~ Q &t any 
macroscopic scale. Let us consider the line Pc = repre- 
sented on Figure 2. In the neutralized model, this path 
crosses a coexistence region (N-L portion). Turning to 
the physical system, the states of this region can not be 
described in terms of phase coexistence, as discussed pre- 
viously : they instead consist in microscopic fluctuations 
constituting a single dishomogeneous phase. In this re- 
strained Pc — subspace, the grand-potential derivatives 
remains continuous {e.g. proton density evolves contin- 
uously) : the physical free energy does not present a 
linear behavior, which means that the phase transition is 
no more of first order. However, the path = is non- 
analytical when it enters the region of pc discontinuity 
(points L and N on Figure 2). The pressure is non- 
analytical at such points, which trace continuous transi- 
tions in the physical system. 

To summarize, as the neutralized-model system enters 
the coexistence region and undergoes a first order phase 
transition, the physical system restrained to the subspace 
Pc = turns into a dishomogeneous phase via a continu- 
ous transition. The only way to keep a first-order phase 
transition in the physical system would be to have no pc 
discontinuity at the phase transition in the neutralized 
system, which could be realized only if the electron dis- 
tribution followed closely the proton distribution leading 
to a complete screening. 

In the neutron-star phenomenology, the mixed phase 
consists in a Coulomb lattice of neutron-rich nuclei at 
zero temperature, or the pasta phase(s) at finite temper- 
ature. Most calculations presented in literature treat the 
electron dynamics at the mean- field level either by con- 
sidering electrons as a uniform background [1, [12, HI], 
or allowing for a slight polarization of the electron gas, 
leading to a screening length of some tens of fm ^ . An 
explicit calculation of the electron screening in the RMF 
coupled to the electric field [1^ confirms this assumption. 
In this condition, the clusters (nuclei or pasta) show a 
net positive charge at their (microscopic) scale, meaning 
that the total charge is an order parameter for the tran- 
sition in the model neutralized system. As a result, the 
physical crust matter cannot be described as a phase co- 
existence : the core-crust transition occurs between uni- 
form matter and a mixed phase presenting microscopic 
dishomogeneities. This transition is necessarily continu- 
ous, contrary to what is almost universally assumed in 
the literature [13, [H [H, Q [ISl. Of course, for such 
dishomogeneous structures to emerge, calculations be- 
yond the mean field are necessary, as will be discussed 
in the next section. In this case the phase diagram can 
be more complicated and present other bifurcations if 
the different pasta structures, here represented by a sin- 
gle coexistence region, correspond to different phases as 
suggested by molecular-dynamics calculations [30|. The 
limited increase of the static form factor to neutrino scat- 
tering observed in the numerical calculations of ref . [l^l is 
consistent with the thermodynamic arguments developed 
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in this section. 



E. Quenching of critical behavior in stellar matter 
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FIG. 4: (Color online) Finite-size scaling for the neutralized 
Ising model with long-range interactions. Upper part : limit- 
ing temperature as a function of the lattice linear size giving 
the exponent v. Medium part : scaling of the order parameter 
giving the consistence between /3 and v. Lower part : scaling 
of the susceptibility giving the consistence between 7 and v. 
The right and left parts of the figure give two different ways 
of fixing the critical exponents. None of them gives a coherent 
scaling of all the quantities. 



Turning now to the mixed-phase phenomenology of 
high-temperature supernova cores, it is interesting to re- 
mark that long-range Coulomb interactions also strongly 
affect critical behaviors. Indeed, the Coulomb-energy 
density can be expressed as a function of charge fluc- 
tuations as 



a 



drdr ~ lixa / ac.[T)rAr , (14) 



where cr(r, f^) = {5pc{r)5pc{r')) is the charge-density 
fluctuation, and translational and rotational invariance 
imposes acir^r') = ac{\r — r'\). If the electron screening 
is a small effect even at high temperature, i.e. if pc keeps 
being an order parameter up to the critical point, then 
<Jc{r) is expected to scale as 



where d is the space dimensionality and 77 is a critical ex- 
ponent which turns out to be close to zero in most phys- 
ical systems (77 = 0.017 for the liquid-gas universality 
class) . A critical point would then correspond to a diver- 
gent Coulomb energy. Hence, only three behaviors are 
allowed. Either the critical temperature goes to infinity, 
or the order parameter becomes independent of the total 
charge or the phase transition in the neutralized system 
ends at a first-order point. The corresponding second- 
order point in the physical system cannot then be a crit- 
ical point, and critical opalescence is suppressed. Con- 
cerning supernova cores, the presence of critical opales- 
cence may have important consequences on the opacity 
to neutrino scattering [l^ , and neutrino transport is cru- 
cial for simulation studies of supernova explosions [1, @| ■ 
According to this qualitative argument, we do not expect 
a strong increase of matter opacity to neutrino scatter- 
ing. To have some quantitative estimation, it is clear 
that calculations of stellar phase structure going beyond 
the mean-field approximation are necessary. 

Our recent finite-size-scaling calculations within a Lat- 
tice Gas model confirm the suppression of critical be- 
havior in stellar matter, as we now show[31]. In this 
schematic but exactly solvable model, each site of a three- 
dimensional lattice oi Q. — sites is characterized by 
an occupation number = 0, 1 and an effective charge 
Qi = n-i — X]jLi"-j/^- The effective charge represents 
the proton distribution screened by a uniform electron 
background. 

The schematic Hamiltonian H = Hm + He with 



Hn 
He 



2^ n^nj 

<i]> 



K gig J 

2 ^ r 



- niHijCij 



(16) 

(17) 



(Jc{f) cx r 



2-d- 



(15) 



is introduced to study the interplay of nuclear- like {Hjq) 
and Coulomb-fike forces {He)- X]<ij> ^ s^'^ extended 
over closest neighbors, and r^j is the distance between 
sites i and j . The short-range and long-range interactions 
are characterized by the coupling constants e and k = 
a.p^'^x^ respectively, where po is the nuclear saturation 
density, and x is the proton fraction. 

To accelerate thermodynamic convergence, the finite 
lattice is repeated in all three directions of space a large 
number iV/j of times. The Coulomb interaction with 
the different replicas of each site is analytically calcu- 
lated and shown to be equivalent to a renormalization 
of the long-range couplings Cij [3l| . The phase diagram 
of the model is evaluated with standard Metropolis tech- 
niques 31] and shown to contain for all simulated lat- 
tice sizes L a coexistence region terminating at a lim- 
iting temperature Tiim.{L). The quenching of critical- 
ity can be formally verified in terms of critical expo- 
nents. If the asymptotic value of the limiting tempera- 
ture Tiim = limi^oo Tiim{L) did correspond to a critical 
point, any generic thermodynamic observable Y charac- 
terized by the critical exponent <i> should fulfil close to 



Y{L,t)^f{L/m)\tf. (18) 

with t = T/Tiini — 1. Scale independence close to the crit- 
ical point imposes for the scaling function /(s) oc s^*/'', 
where v rules the divergence of the correlation length, 
^ oc t~'' . Specifying eg. p^ to the behavior of the limit- 
ing temperature gives 

\fum{L)~Tu^\^L-^'''. (19) 

The scaling of the order parameter Ap — {n)^ ~ {11)0 
and of the susceptibility x — '^ij{5ni5n.j) /T results 

L^^/^Ap cx (ii/^il)'' r<f',„„(L),L>e (20) 

L-'I'^X ^ (Li/^tl)"'' T>fH,„(L),L»e (21) 

while both expressions should become independent of L 
and T as soon as the critical point is approached for each 
given finite size, L ^. For T < TumiL), x contains 
also jumps between the low-density solution {n)c and the 
high-density one {n)^, and should obey first-order scaling 
as t^'^^x oc L'^. Introducing the hyperscaling relation 
d = (7 + 2/3)/t/, we get 

/ , \ 2/3 

L-^/^'x oc [L'/'ltlj T < TumiL),L:$> $ (22) 

Fig. 4 illustrates that finite-size scaling is violated for 
the neutralized long-range Ising model at the approach 
of the limiting temperature: it is not possible to find 
a thermodynamically consistent set of critical exponents 
fulfilling at the same time eqs. (|19l20l21l22p . If we fix 
(3 and i' from eqs. p9l20p (panels a and b in Fig. 4), the 
value of 7 is fixed from the hyperscaling relation. The 
calculations with different lattice sizes do not collapse on 
a single curve in the representation of eas. (|2T|) (Fig. 4c), 
even if the first order scaling eas. f^ is nicely respected 
up to Tiim{L). If conversely we fix 7 and v such as to 
fulfil eg. ipTI) (Fig.4d), this does not affect the behavior 
of the order parameter ea. ((20| (Fig.4e), but the scaling 
of the limiting temperature eq. p^ is violated (Fig.4f). 

The loss of critical behaviour has been already ob- 
served in Ising models with long-range frustrating inter- 
actions, where the coexistence region was seen to end at a 
first-order point . This effect was absent in the mean- 
field approximation and was attributed to fluctuations. 

The fact that no enhancement is observed in the 
static form factor of clusterized matter in the molecular- 
dynamics simulations by Horowitz and coworkers [20| 
may also be an indication of this fiuctuation quenching. 

F. Extension of the pasta phases 

Another remarkable effect that should manifest be- 
yond the mean-field approximation is an extension of the 
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Gas 0.5 Liquid 0.1 



FIG. 5: (Color online) Metropolis calculations of the density 
distributions of the neutralized Ising model with long range 
interactions at different temperatures, for a cubic lattice of 
linear size L = 10 and with a proton fraction Pp/ Pb ~ 0.3 
(cNM). The extension of the coexistence zone, corrected for 
finite size effects through finite size scaling, is also shown. 
The lines labelled nNM give the coexistence zone of neutral 
matter for different lattice sizes, and the extrapolation to the 
thermodynamic limit. 



mixed-phase pasta region in stellar matter with respect to 
standard uncharged nuclear matter. Even if at the end- 
ing point of the coexistence region the correlation length 
may not diverge as discussed in the previous section, we 
can still expect that it will be characterized by large cor- 
related structures. As we infer from eq. pTj) and eq.((H|). 
in the neutralized system the Coulomb energy will be 
minimized in homogeneous partitions corresponding to 
pure phases, and it will be maximal in clusterized par- 
titions. This increase of the energy difference between 
the pure phases and the mixed phase results in an exten- 
sion of the coexistence region, i.e. an extension of the 
dishomogeneous phase corresponding to the mixed par- 
titions situated in between the pure phases for each tem- 
perature. Consequently, the limiting temperature for the 
existence of such mixed phase is increased with respect 
to the uncharged system. This effect has been recently 
observed in the neutralized long-r ang e Lattice-gas model 
presented in the previous section [3l|, as shown in Fig. 5. 
The phase diagram of this model (cNM line in Fig. 5) is 
compared to the one obtained when the Coulomb inter- 
action is neglected, k = in eq. (fT7|) (nNM lines in Fig. 5). 
For both models the phase diagram is extracted from the 
Metropolis simulated density distribution in the grand- 
canonical ensemble, also shown in the figure for different 
chosen temperatures. 

As we have argued in Section IIIDl at the mean- 
field level the two models only differ for the presence 
of the electron pressure in cNM, which is not added in 
Fig. 5 and would transform the coexistence curve into a 
second-order transition line. For a given value of pb, the 
charged-particle density pp = pe = xpb is fixed and so 
is the electron pressure (see eq.®), meaning that the 
microstates explored by the calculation with (cNM) or 
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without (nNM) Coulomb in the canonical ensemble are 
the same at the level of mean field. Fig. 5 shows that this 
is not true once the phase diagram is numerically cal- 
culated from the exact Hamiltonian without any mean- 
field approximation. The limiting temperature decreases 
with the lattice size for nNM fulfilling finite-size scaling 
with standard 3d-lsing exponents [31], while it is almost 
completely independent of L when Coulomb is accounted 
for, reflecting the finite correlation lenght at the limiting 
point and criticality quenching discussed in the previ- 
ous section. More important, the limiting temperature 
is increased for increasing strenght of the Coulomb field 
(about 6% for the proton fraction a; = 1/3 considered in 
Fig. 5). This expansion of the coexistence region is due 
to the energy cost of fragmented inhomogeneous config- 
urations at p/pq ~ 0.5 respect to the the largely uni- 
form proton-charge distribution in pure phases, whose 
Coulomb energy is almost exactly compensated by the 
electron background. This effect is entirely due to the 
screening effect of the electrons, under the assumption 
that they can be approximated by a fixed background. 

As such, the coexistence-region expansion is specific 
of the stellar problem, while other physical systems sub- 
ject to Coulomb frustration exhibit the opposite behav- 
ior. This is notably the case of frustrated Ising fer- 
romagnets, as well as of finite atomic nuclei. In such 
cases the Coulomb repulsion is known to reduce the lim- 
iting temperature [H, HI, [H, [H, [s^. This reduction 
is also a usual expectation in the astrophysical context 
@, S i, El [H, [11 [111- An amplification of the mixed- 
phase region at zero temperature has already been re- 
ported in ref.jl^. In this work a consistent treatment 
of the proton-density rearrangement effect due to the 
Coulomb field under the constraint of strict charge neu- 
trality in the Wigner Seitz cell is seen to de-favour, at 
each given total density, the formation of large charged 



structures (rods respect to droplets, bubbles respect to 
tubes). The net effect is a density amplification of the 
mixed phase. Our results imply that this effect should 
persist at finite temperature. An increase of the max- 
imum temperature for the coexistence region might be 
expected, and the mixed-phase phenomenology may be 
relevant for the proto-neutron-star structure in a wider 
temperature range than usually expected [T6| . 



III. CONCLUSION 

In this paper we have shown that long-range Coulomb 
forces in dense stellar matter require a non-trivial discus- 
sion on the definition of the thermodynamic limit. When 
the constraint of total charge is not strictly zero, the di- 
vergence of the Coulomb energy implies that the chemical 
potential associated with the total charge loses its ther- 
modynamic meaning. Hence, the grand potential, i.e. 
the system pressure, and its derivatives, i.e. the parti- 
cle densities, result independent of the charge chemical 
potential. The suppression of a thermodynamic degree 
of freedom strongly affects the thermodynamics of the 
charged system and the phase-transition phenomenology. 
In general, the Coulomb field modifies the non-analytical 
properties of the partition sums and changes the order 
of the transitions: the first-order core-crust transition 
obtained when the coulomb energy is disregarded, turns 
into a continuous transition from the homogeneous core 
phase to a mixed phase (the so called pasta phases), when 
the coulomb energy is accounted for. A critical point 
can be found at finite temperature if and only if the net 
charge does not present fluctuations. Only in this case, 
critical opalescence in stellar matter can be compatible 
with the long-range Coulomb force. 
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